Optoelectric spin injection in semiconductor heterostructures without ferromagnet 
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We have shown that electron spin density can be gener- 
ated by a dc current flowing across a pn junction with an 
embedded asymmetric quantum well. Spin polarization is 
created in the quantum well by radiative electron-hole recom- 
bination when the conduction electron momentum distribu- 
tion is shifted with respect to the momentum distribution of 
holes in the spin split valence subbands. Spin current appears 
when the spin polarization is injected from the quantum well 
into the n-doped region of the pn junction. The accompanied 
emission of circularly polarized light from the quantum well 
can serve as a spin polarization detector. 
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One of the most important problems in spintronics is 
the efficient injection of spin currents into semiconductor 
structures. A possible way is to inject spin polarization 
from a ferromagnet by passing a dc current across the in- 
terface 1^. If the ferromagnet is metallic, the efficiency of 
such injection is controversial ||^, and the observed weak 
spin current has been attributed to the large mismatch 
of spin diffusion constants between the adjacent semicon- 
ductor and metal Spin injection can be enhanced by 
using an appropriate interface tunneling barrier Q. On 
the other hand, rather high degree of spin current po- 
larization has been detected if the spin is injected from 
magnetic semiconductors although its high efficiency 
is restricted to low temperatures. 

In this Letter we propose a new method to use a dc 
current to inject spin polarization, not from any ferro- 
magnetic material, but from a quantum well (QW) em- 
bedded into a pn junction. The spin polarized current is 
generated during the radiative electron-hole recombina- 
tion in the QW, accompanied by the emission of circu- 
larly polarized light. Our new injection mechanism can 
be easily explained. It is well known that in a QW which 
is asymmetric along the growth direction, at each finite 
wave vector k the spin degeneracy of hole subbands is 
removed. The corresponding splitting of hole energies 
increases with k and can reach a quite high value. For 
example, in a p inversion layer of a GaAs-AlGaAs hetero- 
junction, the splitting of the topmost heavy- hole subband 
was calculated ^] to be about 5 meV at A:=2-10® cm~^. 
Each of the split states is a linear combination of four an- 
gular momentum eigenstates which are specified by the 
z-component J2=±3/2 and ±1/2. The resulting mean 
spins of hole states are parallel to the QW interfaces. 
For a given k the spins in each pair of spin-split sub- 
bands have opposite directions, and in a given spin-split 
subband the mean spin at k is equal and opposite to that 



at -k. The spin orientations in the topmost heavy-hole 
subband are shown in Fig. 1, where the Fermi energy /ig 
(or /i/i) of the quasiequihbrium degenerate electron gas 
(or hole gas) is indicated. When a cr-spin electron in the 
k state of the conduction band makes a radiative tran- 
sition to the topmost heavy-hole subband, as marked by 
the downward vertical arrow line in Fig. 1, the probabil- 
ity of this process depends on a because the population 
of hole at the k state in the split (-l-)-subband may be dif- 
ferent from that in the split (-)-subband. Hence, at k the 
conduction band electron spin will be polarized along a 
direction if an electron with ("(T)-spin has higher recom- 
bination rate. However, since the hole spin orientation at 
-k is reversed with respect to that at k, the spin polar- 
ization of the conduction band electron at the ~k state 
will be along -a direction. Consequently, if the momen- 
tum distributions of both the electron gas and the hole 
gas are isotropic, there will be no net spin polarization 
of the conduction band electrons. 




FIG. 1. Schematic illustration of the energy levels and 
the orientations of mean angular momenta in the topmost 
heavy-hole subband in an asymmetric quantum well. 

On the other hand, if the momentum distributions of 
the electron gas and/or the hole gas are anisotropic, the 
generation of spin polarization becomes possible. We will 
explain the generating process with the help of Fig. I, 
where the anisotropic momentum distribution is indi- 
cated by a shift 5k of the quasiequihbrium momentum 
distribution of electrons with respect to that of holes. 
The state 3 in the (+)-subband and the state 4 in the 
(~)-subband are equally occupied by holes with spins in 
opposite directions. Hence, the total probability of re- 
combination of a conduction electron with the two holes 
in the states 3 and 4 is independent of the spin of the 
electron. However, a hole appears in the state 1 in the 
(+)-subband but not in the state 2 in the (-)-subband. 
Hence, the conduction electrons in the shaded region can 
recombine only with the holes in the (+)-subband, and 



1 



the corresponding recombination probability is spin de- 
pendent. Such processes create a nonequiUbrium electron 
spin polarization in the QW. When this polarization dif- 
fuses out of the QW, a spin current is then generated 
in the n-doped region. One way to realize the situation 
shown in Fig. 1 is to apply an electric field E parallel 
to interfaces. Then, both electrons and holes gain their 
respective drift velocities. Let v be the relative drift ve- 
locity of electrons with respect to drift velocity of holes. 
If we ignore the drift of the low mobile holes, the resul- 
tant band positions are illustrated in Fig. 1. The shift Sk 
is simply m*^ / h=m* fj,E / h, where ^ is the electron mo- 
bility. The suggested device for spin current injector, a 
QW embedded into a pn junction, is illustrated in Fig. 2. 
A bias voltage Vt is applied along the z axis, and the 
transverse voltage source Vt creates an electric field E 
along X axis. The potential profile along z axis and the 
spin current are also shown schematically in Fig. 2. 



tion vector e^, to the lowest order, the quantum proba- 
bility amplitude of this process can be easily derived as 
CiM^^Jk, q)-fC„iM^'|'^(k, q), where z/=± represents 

the (±) split valence subbands, and 

M,^;„-(k,q) = 
pA , 

' dz^;^%^j^{z)xw{z). (1) 
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In the above equation, Pa-.,j^ = {Jz\p\o') is the interband 
matrix element of the electron momentum operator p, 
and Aq=y^2Trc/nq. 

Knowing the recombination probability of the conduc- 
tion band electrons with a specific spin, the spin genera- 
tion during the recombination process can be calculated. 
The resultant spin density in the conduction band is de- 
rived as 



(2) 
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FIG. 2. Suggested device for spin current injector. Upper 
plot is the corresponding potential profile. 

At low temperatures, for degenerate electron and hole 
gases, besides the band to band electron transitions, ex- 
citonic recombination processes must be considered. To 
avoid the complicated analysis which will not change the 
essential physics, in this Letter we will neglect the exci- 
tonic effect on the spin generation. 

For the valence band we define \Jz) as the zone cen- 
ter Bloch state and « j (•^) ^^'^ n-th confined state in 
the QW associated to the two-dimensional wave vector 
k and the hole spin projection J^. Then, the wave func- 
tions '^nk(^) °^ ri-th (±) split valence subbands can 
be expressed as the sum of V'k n ^^^^ '^^-—^1 

and ±i. For electrons in the lowest conduction sub- 
band, we can similarly define \a) and Xki^)- If we ne- 
glect the small spin-orbit splitting, the electron wave 
functions 'l'k(-z) can be written as a linear combination 
of two degenerate states CiX]i{z)\^)+C_iX]iiz)\ — ^), 
where C_(_i are normalized amplitudes. When an elec- 
tron and a hole recombine from their respective quantum 
states to emit a photon of wave vector q and polariza- 



where 2s^r^ are Pauli matrices. The matrix elements 
Gaa' have the form 

G„„,=A M,V;(k,q)M,^r(k,q)/(k)/:(k), (3) 

k,q, A,n,i/ 

where A is a numerical factor, and /(k) [/^(k)] is the 
momentum distribution functions of electrons (or holes). 
Within the perturbation theory, /(k) and /,'^(k) are as- 
sumed to be spin independent. A relevant dimensionless 
parameter which measures the efficiency of spin genera- 
tion is P=Gs/G, where G=Gi i -l-G_ i _ i is the number 

2 ' 2 _ 2 ' 2 _ 

of radiative recombinations per the time unit . 

As mentioned earlier, we will consider higher temper- 
ature such that the excitonic effect can be neglected. 
Hence, each momentum distribution function is the 
sum of the equilibrium Boltzmann distribution and the 
nonequilibrium correction due to the electric field E. 
Since the drift of holes can be ignored, we have 



/„^(k)=exp 
/(k) = exp 
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where e(fc) is the lowest conduction subband, and ^nik) 
(±) split valence subbands. The zero reference energy is 
set at the bottom of the lowest conduction subband. 

We need the valence subband wave functions ip^^^ j {z) 
for calculating of M^'^(k). Let the growth direction z 
be along the [001] crystal. The wave functions can be 
obtained by applying a proper unitary transformation 
1^ which block diagonalizes the Luttinger Hamiltonian. 
They can be expressed in the general form as 

V'^; _3 w = ±'/'L,iW = ±* exp(z</.k)e;^j^)/V2, 

= ±V'^_„_.(z) = ±* exp(zr;k)e^^,(z)/y2, (5) 
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where the real functions ^^/^{z) and ^^i{z) represent the 
partial amplitudes of heavy and light holes in the (+)- 
and the (-)-state. The phases are 0k=3?7k=-3<pk/2 with 
ipii=cos~^ {kx / k) . To obtain (||) we have neglected the 
band warping. For convenience we set the x-axis along 
the drift velocity direction. 

With all above equations (|l|)-(|^) we are ready to cal- 
culate the efficiency of spin generation P=Gs/G. In 
terms of the overlap integrals b^f^=J dz^^f^{z)xi!iiz) and 
bni=I dz^^i{z)x^{z), we define i?^,„,=(&^,)V3, i?±„,=~ 
b^ib^i^/V^, and R^^^=0. Then, P is derived as 



, k,n,z-'— ± 



(6) 



i k,n,z^— ± 



where J^fe^„_y=exp { — [ej'j(fc) -f e{k)]/kBT}. As expected, 
only the anisotropic part (v-k term) of the electron mo- 
mentum distribution function /(k) has contributed to 
the generated polarization. 

Based on the above formula, P can be investigated 
numerically in a broad range of parameters. However, 
qualitative results can be derived from (|^) analytically 
for /cc?<Cl, where d is the typical length of hole confine- 
ment in the QW. We notice that in (||), at a temperature 
T, a characteristic wave vector fc^ can be defined as the 
mean value of k with respect to the Boltzmann factor 



Then, an order of magnitude evaluation of P 



can be obtained at kTd<$^l. Our calculation indicates 
that in the region fcc?<Cl, the energy split A„(A:)=e+(A:)- 
e~(A;) is much less than the separation of two adjacent 
hole subbands, and as well as b^f^2±b~^. Hence, 

R;;;^., and the term {'R;^^Tk,n,+-'R-nk^k,n,-) to be 



summed up in the numerator of (|6|) becomes proportional 
to 1-exp [—An{k)/kBT]. With such simplification we can 
analyze the physical processes which contribute to P. Let 
us consider the contribution of the lowest hole subband 
with n=l. If this is a heavy hole subband, then the small 
admixture of light hole states to this subband allows us 
to evaluate bu to obtain 6i/cx(fcc?)^6ift. As a result, we 
have Ry^ik^Rx.ik, and the generated spin polarization is 
oriented perpendicular to the drift direction of electrons. 
On the other hand, if the n=l subband corresponds to 
the light hole, similar analysis gives bih(x(kd)^bii and 
Rx,ik^Ry,ik- Consequently, the so generated spin polar- 
ization is oriented parallel to the electron drift direction. 

In the range of temperatures where Ai{kT) /ksT^^l^ 
we derive from (|^) 
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hvkx Ai(fcT) 
ksT ksT 



{krdf 



(7) 



where (3—2 (or /3=0) if the lowest subband is of heavy 
(or light) hole type. Hence, in the region of small charac- 
teristic wave vectors, the optoelectric generation of spin 



polarization is more efficient if the lowest hole subband 
is of the light hole type. We should mention that in III-V 
semiconductor QW, the lowest hole subband is of light 
hole type if the QW is sufficiently strained. One example 
is InA^GaAs QW. 

In (0) the factor hvkT/kBT is due to the anisotropic 
momentum distribution of the electrons. The other fac- 
tor Ai {kT){kTd)^ /ksT is originated from the ratio of the 
two summations in (^, and increases with T through kx 
and Ai(A:t). For example, in GaAs at room tempera- 
ture, A;T=i2 • 10^ cm-\ and then fcrrf^l for d=100 A. 
Under this situation, Ai(fcT) becomes comparable to the 
quantization energy of hole subbands and (^ is no 
longer valid. However, a simple scaling analysis shows 
that around kxd^l, Eq.(0) can still be used to evaluate 
V with the factor {kxd)^ replaced by a numerical factor 
of the order of unity. 

We have shown that spin polarization of conduction 
electrons can be generated in a quantum well due to the 
radiative electron-hole recombination. Let \e\I be the 
electric current across the pn junction, and ri=Tnr/ [Tnr + 
Tr) the luminescence quantum efficiency, where (or 
Tnr) is the radiative (or nonradiative) electron-hole re- 
combination time. Then, the number of radiative recom- 
binations per time unit is G=/r/. The so generated spin 
polarization in the conduction band can either diffuse out 
of the QW into the n-doped region of the pn junction, 
or relax within the QW with a spin relaxation time Tsw ■ 
In III-V semiconductor QWs the dominating process is 
the the D'yakonov-Perel' spin relaxation although the 
Bir-Pikus mechanism j|] can be efficient due to presence 
of a large number of holes. In the steady state, the spin 
current which flows out of the QW is given as 



Is = VIv - S/t - S/t, 



(8) 



where S is the spin polarization, and 1/T=l/Tr+1/T„r- 
We choose the spin quantization axis parallel to the di- 
rection of the vector P defined in (^) . Then, S is simply 
(ni-n_i)/2, where n„ is sheet electron density in the 
QW with cr-spin. 

To relate S to /, we need to specify the transport pro- 
cess between the QW and the n-doped region, which in- 
volves the spin diffusion and relaxation in the n-doped 
semiconductor. To avoid such complication, we will con- 
sider a simple model of thermionic transport over the 
barrier, which is suitable for not too low temperatures. 
In this case the cr-spin current v is determined by the 
balance of Richardson currents emitted from both sides 
of the potential barrier. Such emission currents depend 
on the chemical potentials pL^w in the QW and h^ 
the n-doped bulk semiconductor. In the linear response 
regime, \^icrw-IJ^(jb\-^kBT , and we have 
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(9) 



where A* is the Richardson constant, and the the barrier 
height U is indicated in Fig. 2. With this expression, the 
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spin current Is = {ii-i_i)/2 is a function of the chemical 
potentials. 

In the n-doped bulk semiconductor we assume diffu- 
sive motion of electrons with a diffusion constant D and a 
spin relaxation time Ts&. We should mention that within 
the framework of linear transport theory the spin current 
in a bulk material can be driven only by a spin density 
gradient, and the corresponding characteristic length of 
the spin density variation has the form Lg~^ Drsb- In 
an n-doped bulk semiconductors Tgh can be very long , 
and a large Ls will reduce the efhciency of spin injec- 
tion. This is the same problem encountered in the study 
of spin injection from ferromagnetic metals into semi- 
conductors 1^, which can be overcome with a proper 
choice of the barrier height U . If U is so chosen that 
(Ls/Z) exp [(/io-fe — eJ7)/A;BT]^l, where I is the electron 
mean free path, the injected spin current is no longer 
limited by the low spin relaxation rate of the bulk semi- 
conductor. For a moderately doped III-V semiconduc- 
tors, if we take from [||J|] rsb=100 ps and the electron 
mobility 5-10'*cm^/Vsec, Ls//~10 at room temperature. 
Therefore, the above inequality can be easily satisfied 
even if the barrier height is rather low. In this situation 
we obtain 
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ksT Tsw ' 



(10) 



where A/i=(^ i i ^^-^i i _f,-M-i,f,)/2- 

Although ^yujkBT is assumed to be small in linear 
response theory, the ratio t/tsw can be large, and so 
the second term in the square brackets of ( |lO[ ) can not 
be neglected. For example, in bulk III-V semiconduc- 
tors with carrier density in the range between 10^^ cm"'^ 
and 10^^ cm"'^, at room temperatures the ratio t/t^w 
is about ten |^,^, assuming t in the range of nanosec- 
ond. The similar ratio was also found in QWs jl^. In 
our system, as illustrated in Fig 2, the electron density 
and the hole density in the asymmetric QW are spatially 
shifted with respect to each other. The electron-hole re- 
combination becomes spatially indirect, and therefore r 
increases. On the other hand, Tsw increases if the spin 
relaxation is dominated by the Bir-Pikus mechanism, but 
remains almost the same if the D'yakonov-Perel' mecha- 
nism dominates. Consequently, the actual value of t/t^w 
varies with the system to be studied. From (|l^), we see 
that the upper limit of the spin polarization of the in- 
jected current is Ig/I—Vr]. 

Since the spin polarization in the QW is created by 
the emission of circularly polarized light, the spin gen- 
eration can be detected by measuring the polarization 
of emitted photons. Circular polarization of a photon 
is represented by the imaginary off-diagonal elements of 
the photon polarization matrix p. This Hermitian matrix 
can be calculated in a way similar to the calculation of 
spin polarization. It can be shown that the imaginary ele- 
ments are related to P as pxzC<.iPy and pyzOH-iPx- Hence, 
most of the circularly polarized photons are emitted with 



their wave vectors parallel to P. 

If we introduce a ferromagnetic layer to the n-doped 
part of our sample, the optoelectric injection of spins 
can be investigated by measuring the resistance of the 
system. The separation of the ferromagnetic layer and 
the QW must be less than Ls- In this case, the measured 
resistance depends on relative spin polarizations in the 
paramagnetic and the ferromagnetic materials |^,|^ , and 
can be varied by changing the direction of the applied 
electric field E. 

We close this Letter with one remark. Within the 
framework of linear response theory used in our analy- 
sis, the ratio hvkT/ksT must be small. With a stronger 
electric field to increase the electron drift velocity v, we 
conjecture an enhancement of the spin generation under 
an anisotropic momentum distribution of hot electrons. 
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